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Universal consequences of the presence of excluded volume interactions in dilute 

polymer solutions undergoing shear flow 

K. Satheesh KumaiQ and J. Ravi Prakash 
Department of Chemical Engineering, Monash University, Clayton, Victoria-3168, Australi^\ 

(Dated: February 2, 2008) 

The role of solvent quality in determining the universal material properties of dilute polymer 
solutions undergoing steady simple shear flow is examined. A bead-spring chain representation of 
the polymer molecule is used, and the influence of solvent molecules on polymer conformations is 
modelled by a narrow Gaussian excluded volume potential that acts pair-wise between the beads 
of the chain. Brownian dynamics simulations data, acquired for chains of finite length, and extrap- 
olated to the limit of infinite chain length, are shown to be model independent. This feature of 
the narrow Gaussian potential, which leads to results identical to a 5-function repulsive potential, 
enables the prediction of both universal crossover scaling functions and asymptotic behavior in the 
excluded volume limit. Universal viscometric functions, obtained by this procedure, are found to 
exhibit increased shear thinning with increasing solvent quality. In the excluded volume limit, they 
are found to obey power law scaling with the characteristic shear rate f3, in close agreement with 
previously obtained renormalization group results. The presence of excluded volume interactions is 
also shown to lead to a weakening of the alignment of the polymer chain with the flow direction. 



I. INTRODUCTION 

The nature of the thermodynamic interactions be- 
tween the monomers in a polymer chain and the solvent 
molecules in their neighborhood, determines the ensem- 
ble of spatial configurations adopted by the chain, and as 
a result, has a significant influence on all conformation 
dependent properties of polymer solutions. Extensive in- 
vestigations of static polymer solutions has established 
that both the temperature of the solution and the molec- 
ular weight of the dissolved polymer, control the strength 
of these interactions. Indeed their influence can be com- 
bined together in terms of a single variable, the so-called 
solvent quality parameter, z = vo(l — Tg /T)\J~M \ where, 
vq is a polymer-solvent system dependent constant, T is 
the temperature of the solution, Tg is the theta tempera- 
ture, and M is the molecular weight. At large values of z, 
usually termed the excluded volume limit, various proper- 
ties of polymer solutions have been shown to obey power 
law scalings with molecular weight. On the other hand, 
at intermediate values of z (the crossover region), the 
behavior of a host of different polymer-solvent systems 
has been found to be collapsible onto universal master 
plots with a suitable choice of the constant vo i 1 ' 2 '?' 4 

The entire range of observed behavior has been suc- 
cessfully modelled by static theories of polymer solutions 
by mimicking the influence of the solvent molecules on 
the conformations of a polymer coil, with the help of a 
potential that acts pair-wise between the various parts 
of a polymer chain. Typically, by representing the linear 
polymer molecule by a bead-spring chain consisting of N 
beads connected together by (N — I) Hookean springs, 
more or less repulsion by the solvent molecules is ac- 
commodated by varying the magnitude of the pair-wise 
excluded volume interaction parameter, denoted here by 
z* , between the beads^Si^ In particular, by noting that 
z* ~ (1— Tg/T), and N ~ M, a quantitative prediction of 
both the excluded volume limit and the crossover regime 



behavior of a number of equilibrium properties has been 
achieved by analytical theories, with the help of the map- 
ping z — z*^/N. 

The growing recognition that polymer-solvent inter- 
actions also play a key role in determining the behav- 
ior of polymer solutions far from equilibrium is respon- 
sible for the recent appearance of a number of studies 
aimed at obtaining a better understanding of this phe- 
nomenon, 9 ! 10 . 11 ! 12 . 13 ! 14 . 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 Due to the an- 
alytically intractable nature of the non-equilibrium the- 
ory, a majority of these investigations have been largely 
numerical in nature. As in the case of static theories, 
polymer-solvent interactions have been modelled with 
the help of an excluded volume potential, such as the 
Lennard-Jones potential, acting pair-wise between the 
different parts of the polymer chain. While a number 
of qualitative conclusions on the influence of excluded 
volume interactions have been drawn by these studies, 
a methodical program aimed at cataloging the univer- 
sal behavior of polymer solutions in a variety of flow 
fields in terms of the solvent quality parameter z — as has 
been done at equilibrium — has not yet been attempted. 
Nevertheless, some progress in this direction has been 
made in the framework of approximate analytical the- 
ories. In the excluded volume limit, Ottinger and co- 
workersSiiSiS^ have examined, with the help of approxi- 
mate renormalization group methods, the universal be- 
havior of polymer solutions undergoing shear flow. The 
solvent quality crossover scaling behavior of viscometric 
functions was, however, not reported by these authors. 
Recently, Prakasbi 7 ^ 9 . has obtained, with the help of a 
Gaussian approximation, universal material functions in 
shear flow, both in the crossover regime and in the ex- 
cluded volume limit. 

A systematic methodology by which it is possible to 
go beyond approximate results, and obtain exact results 
(albeit within numerical error bars) of the influence of 
solvent quality, in terms of the parameter z, on both 
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equilibrium and non-equilibrium properties, has been 
newly introduced by Kumar and Prakash^4 Basically, 
they showed that model independent predictions can be 
obtained by extrapolating exaci^ Brownian dynamics 
simulations data, acquired for chains of finite length, to 
the limit of infinite chain length. While their method- 
ology is applicable to both regimes of behavior, they re- 
stricted their attention to the prediction of equilibrium 
and linear viscoelastic properties. In this paper, we ap- 
ply their procedure and obtain, for the first time, exact 
predictions of the universal crossover dependence of vis- 
cometric functions on z, at various values of the charac- 
teristic non-dimensional shear rate, f3 = A P 7, where A p 
is a characteristic relaxation time defined by, 



P " N A k B T 



(1) 



Here, &b is Boltzmann's constant, ?7 S is the solvent vis- 
cosity, [77] o is the zero shear rate intrinsic viscosity, and 
JVa is Avogadro's number. The same methodology is also 
used here to calculate the dependence of several proper- 
ties on (3, in the excluded volume limit. This includes, 
for instance, the power-law scaling of viscometric func- 
tions with (3 for large values of j3, and the dependence 
of the orientation of the polymer coil with respect to the 
flow direction, on the magnitude of f3. The important 
question of whether the value of the critical exponent v 
(which governs the scaling of the mean size of a polymer 
chain with molecular weight in the long chain limit), is 
unaltered in shear flow, is also addressed here within this 
framework. 

In addition to excluded volume interactions, it is now 
well established that it is also essential to include sol- 
vent mediated hydrodynamic interactions between dif- 
ferent segments of a polymer chain, in order to obtain an 
accurate description of the universal dynamics of poly- 
mer solutions. While several studies (based on bead- 
spring chain models) of the combined influence of ex- 
cluded volume, hydrodynamic interaction, and even finite 
extensibility effects, on rheological properties have been 
reported so far, they have generally been restricted to 
finite chainsi 16 i 21 i 22 i 25 i 26 However, Zylka and OttingeriS 
have examined the universal consequences of the presence 
of both excluded volume and hydrodynamic interaction 
effects on properties far from equilibrium, with the help 
of approximate renormalization group methods. Several 
important advances in recent years, within the frame- 
work of Brownian dynamics simulations, have now made 
it feasible to treat hydrodynamic interaction effects, free 
of any simplifying approximations, in bead-spring chain 
models with a large number of beads^i 2 ^ In this work, 
however, we confine attention to excluded volume inter- 
actions alone, for two reasons. First, we wish to exam- 
ine, in isolation, the influence of excluded volume inter- 
actions on long chain properties, before considering the 
more complex possibilities that arise in the presence of 
other non-linearities. The preliminary results of Prab- 
hakar and Prakash. 21 who show that in fact the various 



non-linearities appear to be decoupled from each other, 
suggests that there is some value in looking at these ef- 
fects in isolation from each other. Second, and more prac- 
tically, the procedure adopted here of extrapolating finite 
chain simulation data to the long chain limit becomes 
significantly more computationally intensive when both 
excluded volume and hydrodynamic interactions are in- 
cluded. The results of the present paper are therefore 
not yet directly comparable with experiments. Never- 
theless, while the more ambitious task of a complete de- 
scription of polymer solution dynamics is being simulta- 
neously pursued, it is felt that several important conclu- 
sions with regard to the role of solvent quality, can still 
be made within the present approach. 

The plan of the paper is as follows. In scc.[n] we state 
the governing equations, define the various properties cal- 
culated in this work, and summarize the previously intro- 
duced procedure for obtaining exact predictions of prop- 
erties in the crossover regime and in the excluded volume 
limit. In sec. IIII Al asymptotic results describing the de- 
pendence of the normalized viscosity and the normalized 
first normal stress difference on (3 and z, are presented. 
The excluded volume limit dependence on /?, of various 
properties, is discussed in sec. IIII Bl along with a compar- 
ison with results obtained previously with the Gaussian 
approximation. Finally, the principal conclusions of this 
work are summarized in sec. IIVI 



II. 



BEAD-SPRING CHAINS WITH EXCLUDED 
VOLUME INTERACTIONS 

A. Governing equations 



Within the framework of a numerical solution strategy, 
the most appropriate model for the task of describing the 
rheological behavior of dilute polymer solutions, in the 
limit of long chains, is a bead-spring chain model. The 
internal configuration of a bead-spring chain, suspended 
in a Newtonian solvent undergoing homogeneous flow, 
with the beads located at positions r u , v = I , . . . , N with 
respect to an arbitrarily chosen origin, can be specified 
by the N—l bead connector vectors, Qk = r/,+1— rj., k — 
1, . . . , N— 1, connecting beads k and k + 1. For Hookean 
springs, and a pairwise intra-molecular excluded volume 
force, the potential energy of a bead-spring chain is given 
by the expression, 



1 N-l N 

-H J2 Q<-Q< + * E 



E(r v 



(2) 



V, ^L — l 



where, H is the spring constant, and E (r v — r M ) is the 
excluded volume potential between the beads v and fi 
of the chain. The key ingredient, in the procedure for 
obtaining universal predictions in terms of the parameter 
z, is the use of the narrow Gaussian potential to represent 
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excluded volume interactions, 



components satisfy, 



E(r v 



k B T exp • 



H 



2k B T d* 



(3) 



(W* j (t*))=0 
(W^(t*)W^ k {t*')) = min(f,f')^(S, 



(6) 



where, r vll = r v — r M , is the vector between beads /i 
and v, the nondimensional parameter z* , as mentioned 
previously, measures the strength of the excluded volume 
interaction, and d* is a measure of the range of excluded 
volume interaction. Prakash and coworkers have shown, 
in their examination of universal linear viscoelastic and 
steady shear flow properties in the context of the Gaus- 
sian approximation^*^ and in their calculation of uni- 
versal equilibrium and linear viscoelastic properties with 
Brownian dynamics simulations^ that the usefulness of 
the narrow Gaussian potential stems from the fact that, 
as N — > oo, (i) the parameters z* and N combine to- 
gether to form the single variable, z — z*\~N : and (ii) 
the choice of d* is inconsequential because the parameter 
d* always appears in the theory as the re-scaled variable 
d* /\J~N. These results form the basis of our exploration, 
in this work, of the universal behavior of bead-spring 
chains subject to shear flow. 

As will be discussed in greater detail shortly in sec. Ill Bl 
below, all properties of interest in the present work can 
be calculated from appropriately defined configurational 
averages. Unfortunately, it is currently infeasible to solve 
the nonlinear diffusion equation that governs the time 
evolution of the configurational distribution function in 
the presence of excluded volume interactions. However, 
configurational averages can be numerically evaluated di- 
rectly from ensembles of polymer configurations. Basi- 
cally, stochastic trajectories that describe the temporal 
evolution of an ensemble of chain configurations can be 
generated by solving the stochastic differential equation 
that governs the bead connector vectors, with the help of 
Brownian dynamics simulations^* In terms of the non- 
dimensional variables, 



Qt 



i 



t 



(4) 



where, I = ^k B T/H is a length scale that is propor- 
tional to the root mean square extension of a single spring 
in the Rouse model, Ah = £/4if is a time constant (with 
£ representing the bead friction coefficient), and n(t) is 
the traceless transpose of the velocity-gradient tensor for 
the homogeneous flow field (which can be a function of 
time but is independent of position) , the stochastic differ- 
ential equation that governs the trajectories of the bead 
connector vectors in the presence of excluded volume in- 
teractions, is given by the expression^ 



dQ* 



N-l 



k=l 



dQt 



[7 N 



(5) 

Here, W* is a Wiener process, whose 37V— dimensional 



for j. k — 1,2,3 and is, fi — 1, 2, 3, . . . , N. The quantity 
Bkv is an (N — 1) x N matrix defined by, B^v — $k+i, v — 
<5fc„, with <5fc„ denoting the Kronecker delta, and Ajk is 
the Rouse matrix, 



N 
v=l 



2 for \j - k\ = 0, 
-1 for|j-fc| = l, (7) 
otherwise 



Before discussing the two different simulation schemes 
used here for solving Eq. JSJ to obtain predictions of 
the universal crossover and excluded volume limit behav- 
iors, it is appropriate to first introduce all the properties, 
and their defining equations, that are calculated in the 
present work. 



B. Material functions 

We confine our attention in this work to steady simple 
shear flows, which are specified by the following matrix 
representation for the tensor k in a laboratory-fixed co- 
ordinate system, 



K 



'0 1 0^ 
7 I 
,0 0, 



(8) 



where, 7 is the constant shear rate. 

The principle material functions of interest here are 
the polymer contribution to the viscosity, r/p, and the 
first normal stress difference, vPi. Since hydrodynamic 
interactions have not been included in the present model, 
the second normal stress difference is identically zero^ 
In terms of the components of the polymer contribution 
to the stress tensor, r^ yl T^ x , etc, i] p , and ^1, can be 
obtained from the following relations^ 



r p = 

xy 



r p = 
yy 



-7 2 *i 



(9) 



Within the framework of polymer kinetic theory, the 
polymer contribution to the stress tensor, t p , can be 
calculated by either the Kramers or the Giesekus expres- 
sions. While the latter expression becomes invalid in the 
presence of hydrodynamic interactions, it has been shown 
to remain valid in the presence of excluded volume inter- 
actions^ The relative magnitude of the variance associ- 
ated with the use of either of these expressions in Brown- 
ian dynamics simulations has been examined previously 
by Prakash^ and it has been shown that the Giesekus 
expression always leads, particularly at low shear rates, 
to a smaller value of variance. For this reason, we have 
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used the Giesekus expression in our work, which at steady 
state has the form, 



N-l 



^ ' G m n { K ' {QmQn) + {QmQn) ' 

m,n— 1 

(10) 

where, C mn is the Kramers matrix. The Kramers matrix 
is the inverse of the Rouse matrix, and is defined by, 



N 

C mn = B vm B vn = mm (to, n) — 



(11) 



Here, the N x (N — 1) matrix B v ^ is defined by, B v ^ = 
k/N — O (k—v), with (k — v) denoting a Heaviside step 
function. As will be discussed in greater detail in sec. Ill CI 
below, we further reduce the variance in our simulations 
by implementing a variance reduction procedure based 
on the method of control variatesi^ 

In the limit of zero shear rate, the ratio U^, n , which is 
defined by the expression, 



9 



(12) 



where, n p is the number density of polymers, has been 
shown to have a universal value, both in theta solu- 
tions?^ and in good solvents^ In the context of a dumb- 
bell model, Prakash and Ottinger— have evaluated its 
dependence on the non-dimensional shear rate XhI- In 
this work, we seek to examine the universal dependence 
of Uy v on the characteristic shear rate /3, in the excluded 
volume limit. 

In addition to rheological properties, it is also of inter- 
est to examine the anisotropy induced in any arbitrary 
tensorial property of the solution, due to the orientation 
and deformation of polymer coils caused by the flow field. 
The precise definition of anisotropy is taken up shortly 
below. Of particular relevance in this work, is the change 
in the degree of anisotropy brought about by the presence 
of excluded volume interactions. 

In order to obtain an estimate of the anisotropy of 
any tensorial quantity er, we begin by noting that in 
simple shear flows er must reflect the symmetry of the 
flow field. In particular, since shear flows are invariant 
when the direction of the z-axis is reversed^ the xz, yz, 
zx and zy components of er must be zero. As a result, 
er must have the following matrix representation in the 
laboratory-fixed coordinate system, 



<7= er 



It is clear from Eq. I|13|l that the z-axis is a principal di- 
rection. It follows that the other two principal directions 
must lie in the xy plane. We denote the principal axes in 
this plane by x and y, with the x axis making an angle \ a 



with the i-axis. Thus by rotating the coordinate system 
by Xa-, the matrix that represents ercan be diagonalized. 
By exploiting the rules for change of tensor components, 
the orientation of the (x, y, z) coordinate system relative 
to the (x, y, z) coordinate system can be calculated to be, 



tan 2xo 



^^xy 



(14) 



With the help of symmetry arguments, one can show that 
the orientation angle \<y = t/4 close to equilibrium. As 
the shear rate increases, x<j is expected to decrease to 
zero, reflecting the increasing alignment of the polymer 
coils with the flow direction. The orientation angle x CT , 
therefore, serves as a measure of the anisotropy of er>i£ 

It is straightforward to show that the orientation angle 
associated with the polymer contribution to the stress 
tensor, Xt, is given by ; 15 i 32 



tan 2\ T 



^ T xy 



(15) 



where the quantity to t , defined by the above expression, 
is frequently referred to as the orientation resistance. 52 
In the Rouse model, m T is a constant, independent of 
the shear rate, with a value 2.5. In a similar manner, the 
orientation angle associated with the radius of gyration 
tensor, G, is given by the expressionfi^S 



tan 2\G 



2G 



xy 



Gy 



where, tog is the corresponding orientation resistance, 
and G xy etc., represent components of the tensor G, de- 
fined by, 



G 



1 



JV 

TV ^ 

i/=i 



(ir v - r c )(r v - r c )) 



(17) 



Here, r c denotes the position of the center of mass. In 
the Rouse model, tog = 1.75. It is worthwhile noting 
that while the orientation of the macromolccule mea- 
sured by flow birefringence coincides with Xt (according 
to the "stress optical" law), xg is the orientation an- 
gle measured by static light scattering^ In this work, 
we show that in the presence of excluded volume inter- 
actions, both to t and tog are not constant, but rather 
are functions of the characteristic shear rate /?, and we 
examine this dependence in the excluded volume limit. 
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C. Simulation schemes 








CTzz) 







Two different simulation schemes have been introduced 
and discussed in detail by Prakash and co-workers in or- 
der to find the universal properties of polymer solutions 
in the crossover region, and in the excluded volume limit, 
respectively^2i2i Here, we summarize them briefly. 
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The universal dependence of various properties on the 
parameter z in the long chain limit is found by a two 
step procedure: (i) Simulations are carried out for in- 
creasing values of N, keeping the value of z (= z*^/N), 
constant, and (ii) The accumulated finite chain data for 
any property is then extrapolated to the limit N — > oo. 
The crossover behavior of the property is obtained by 
repeating steps (i) and (ii) for a number of values of z. 
In the Gaussian approximation, the value of the param- 
eter d* was held constant during the process of accumu- 
lating finite chain data^i^i^ and it was shown that the 
particular choice of value for d* was inconsequential in 
the limit TV — > oo. On the other hand, while obtaining 
the crossover behavior of equilibrium and linear viscoelas- 
tic properties with Brownian dynamics simulations^ the 
parameter d* was chosen such that, 

d* = k(z*) 1/5 (18) 

where, k is an arbitrary constant. Note that, for a fixed 
value of z, since z* — > (or equivalently, T — > Tg), as 
N — > oo, this choice of d* implies that the asymptotic 
limit is reached along trajectories in the (z*,d*) param- 
eter space that converge to the origin. It was shown by 
Kumar and Prakash^i that choosing d* values accord- 
ing to Eq. H18p. permits the use of larger step sizes in 
the numerical integration scheme. The independence of 
the asymptotic results from the particular choice of fc 
(and consequently the trajectory in the (z*,d*) param- 
eter space), was also established in ref. I24J. Since Brow- 
nian dynamics simulations are used in the present work 
to obtain the crossover behavior, d* values will be chosen 
according to Eq. (|TH|) . 

Universal properties in the excluded volume limit (z — > 
oo) are obtained by accumulating finite chain data for 
constant values of z* (which is equivalent to the tem- 
perature T being constant), followed by extrapolation to 
the infinite chain length limit. Since for any constant 
value of z*, z — > oo, as N — ► oo, we expect that prop- 
erties obtained in the excluded volume limit will be in- 
dependent of the particular choice of z*. Prakash and 
co-workersi2i2i have shown that this is indeed the case, 
for the situations examined previously by them. As will 
be discussed in greater detail subsequently, this indepen- 
dence will be used in the present work to indicate the 
adequacy or otherwise of the accumulated finite chain 
data. 

The various configurational averages required to cal- 
culate the properties listed in sec. Ill Ul above were ob- 
tained with a second order predictor-corrector Brown- 
ian dynamics simulation algorithm originally proposed 
by Iniesta and Garcia de la Torres A detailed discus- 
sion of the algorithm, adapted for the problem at hand, 
can be found in ref. Il9i In a typical simulation, roughly 
80,000 trajectories, with Gaussian initial distributions, 
were generated, and the integration was carried out un- 
til a stationary state was attained. Message Passing In- 
terface (MPI) clusters were used to distribute the tra- 
jectories across several processors. Since the present 



algorithm belongs to the class of embarrassingly par- 
allel algorithms, the wall time scales with the number 
of processors used. Typically, three different time steps 
At* = 1.0, 0.6, and 0.5 were used, and the results were 
then extrapolated to zero time step using the subrou- 
tine TEXTRA suggested by Ottinger— . This extrapolation 
procedure enabled the simulations to be carried out at 
relatively large time steps. 



A key ingredient in the simulation strategy adopted 
here, as mentioned perviously, is the use of a variance 
reduction procedure based on the method of control vari- 
ates, in order to reduce the statistical error in the simula- 
tions. This is crucial because the magnitude of the error 
in the finite chain data significantly affects the accuracy 
of the extrapolation to the infinite chain length limit. 
Since a fairly involved discussion of the procedure has 
been given previously in ref. l24l it is not repeated here. 
However, a few salient features are highlighted below. 



The essence of the scheme (and the vital factor that de- 
cides its success or otherwise) is to find a control variable 
whose mean value is known exactly, and whose fluctua- 
tions are correlated with the variable whose variance we 
are interested in reducing. It was found by Kumar and 
PrakashjSi that the Rouse model serves as an excellent 
source of control variables for all the equilibrium and lin- 
ear viscoelastic properties that were evaluated by them. 
In this work, we have again used the Rouse model as the 
source of control variables, and found that as before, this 
leads to a significant reduction in the variance of all prop- 
erties. At first sight, this result is somewhat surprising. 
In the case of equilibrium and linear viscoelastic prop- 
erties, since only equilibrium averages were evaluated, 
both the stochastic differential equation, Eq. [S] and the 
stochastic differential equation in the Rouse model, were 
solved with the term involving k set equal to zero. In 
the present instance, since we are interested in properties 
at finite shear rates, both equations were integrated with 
the k term in place. In spite of including this term in the 
Rouse model, as is well known, viscometric functions pre- 
dicted by it are constant and independent of shear rate. 
This is not the case, however, in the presence of excluded 
volume interactions, where, as we shall see shortly, visco- 
metric functions decrease with an increase in shear rate. 
The success in the variance reduction strategy suggests, 
therefore, that even though the Rouse viscometric func- 
tions are independent of shear rate, the fluctuations in 
the Rouse control variables remain correlated with their 
corresponding variables in the presence of excluded vol- 
ume interactions. The reason for this occurance might 
lie in the fact that with increasing shear rate, the beads 
of the chain are drawn further apart, and the resultant 
weakening of the excluded volume interactions actually 
brings the results of the two models closer together. 
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III. UNIVERSAL PROPERTIES 
A. Crossover behavior 

In this subsection, we are concerned with obtaining 
the crossover behavior of the universal ratios, rj p /r] p ^, 
and \& 1/^1,0 (where, r] p fi and ^1,0 are the zero shear 
rate polymer contributions to the viscosity and first nor- 
mal stress difference, respectively), as a function of the 
characteristic shear rate (3, using the simulation scheme 
described in sec. Ill CI above. A strict pursuance of this 
procedure would require the construction of each of these 
ratios at various values of N, using equilibrium simula- 
tions (as described previously in ref. l24|) to find rj p fi and 
\?i,0) and simulations at finite shear rates (as described 
above) to find r] p and ^1, followed by extrapolation to 
the long chain limit. Unfortunately, this procedure suf- 
fers from a serious problem. Since both the denominator 
and the numerator are found by Brownian dynamics sim- 
ulations, the ratios have relatively large statistical error 
bars, which makes their extrapolation to the N — > 00 
limit highly inaccurate. Fortunately, this difficulty can 
be overcome by adopting the following procedure. 

We first note that. 



limjv-i-oo^p/^) 
lim W ^oo (v P ,o/Vp) 

iV^c * li0 lim W ^oo (*i, /*f) 



lim 



lim 



Vpfi 



(19) 
(20) 



where, r] p and \E'f i are the viscosity and first normal stress 
difference predicted by the Rouse model, respectively^ 



■q p = n p \ H k B T 
*f = 2n p X 2 H k B T 



(N 2 - 1) 



3 

(N 2 



1)(2N 2 + 7) 



45 



(21) 
(22) 



Kumar and Prakash have previously evaluated the ra- 
tios rj Pi o/ri p and ^/i^/^f in the long chain limits For 
the purpose of evaluating rj p /rj p xi in the long chain limit, 
therefore, it suffices to evaluate the ratios, rj p /r) p and 
^i/'i'f, as N — ► 00. The simulation scheme described 
in sec. Ill CI above has, consequently, been used here to 
evaluate these ratios at various values of z. 

Figures Q] clearly indicate the independence of the ex- 
trapolated values of r] p /r] p and , in the N — > 00 
limit, from the trajectories in the (z*,d*) parameter 
space used to obtain them. The parameters z and /? have 
been held constant at a value of unity, while carrying out 
simulations at increasing values of N, for three different 
values of k. In each simulation, z* = z/y/N, while d* is 
found from eq 1181 In order to keep (3 constant in each 
simulation, the non-dimensional shear rate Xh 7 used in 
the simulation, at each value of N, has been calculated 
from the expression, A# 7 = j3 (XH n pkBT/rj p ), which 
can be derived from Eq. ^for a dilute polymer solution. 
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FIG. 1: (a) The ratio of the viscosity to the Rouse model vis- 
cosity, and (b) The ratio of the first normal stress difference to 
the Rouse model first normal stress difference, versus l/\fN 
at z — 1 and j3 = 1, for three different values of k. The sym- 
bols, with error bars, represent the simulation results, while 
the lines are least-square curve fits to the data. 



As mentioned earlier, (rj Pj0 / X^ripkBT) can be evaluated 
by carrying out equilibrium simulations for the same set 
of parameter values as those used in the finite shear rate 
simulations. Extrapolation of the finite chain data has 
been carried out by plotting simulation results versus 
1 / \/N, since Prakash 1 ^ has shown previously that lead- 
ing order corrections to the infinite chain length limit, 
of various material properties, are of order l/yN. Fur- 
thermore, as mentioned earlier, the parameter d* always 
occurs in the theory as the ratio d* /\/N. The asymp- 
totic values obtained in this manner correspond to the 
viscosity and first normal stress difference ratios for a 8- 
function excluded volume potential, at z = 1 and j3 = 1. 
The same procedure was then repeated for various values 
of z and (3. 

The universal dependence of r) p /r) Pj o, and ^i/^i.o on 
j3, at various values of z, obtained by combining the infi- 
nite chain length results for the ratios r\ P lf} p and 
with the infinite chain length zero shear rate ratios evalu- 
ated previously by Kumar and Prakash^* are presented 



7 



- 


* 

$ * 

R 




i 

(a) ; 

- 


- 

- 




* 


- 
- 


- 




8 


* ~ 


1 

T 


Z= 1 .U, K= 1 .U 




* 


r x 


z= 1.0, k=1 .5 






□ 


z=5.0, k=1.5 




W 


- o 


z=5.0, k=2.0 






A 


z=10.0, k=1.0 




1 : 


V 


z=10.0,,k=1.5 




i 




1 




10 



* 



+ 


z= 


1.0, k= 


=1.0 


X 


z= 


1.0, k= 


=1.5 


□ 


z= 


5.0, k= 


=1.5 


o 


z= 


5.0, k= 


=2.0 


A 


z= 


10.0, k= 


=1.0 


V 


z= 


10.0, k= 


=1.5 



□ 

A 
V 



(b) 



10 




FIG. 2: (a) Asymptotic nondimensional viscosity ratio 
(Vp/Vp.o), an d (b) Asymptotic nondimensional first normal 
stress difference ratio (^i/^i^), versus the characteristic 
shear rate (3, for various values of z and k 



in Figs. [3 At each value of z, extrapolated data for 
two different values of k are displayed in order to clearly 
delineate the limit of (3 up to which the current asymp- 
totic results are valid. In essence, beyond some threshold 
value of the characteristic shear rate, say /?', the lack of 
coincidence of data for two different values of k, implies 
that the present finite chain data, accumulated for chains 
with N < 36, is no longer sufficient to obtain an accu- 
rate extrapolation. For f3 > (3\ the flow field begins to 
probe model dependent length scales, and scale invari- 
ance, which is responsible for the observed model inde- 
pendence, no longer exists. The results in Figs. suggest 
that /3' for V&i/^i.o is smaller than that for r] p /r]pfi, an d 
in both cases, it decreases with increasing z. 

It is clear from Figs. |2] that the presence of excluded 
volume interactions leads to shear thinning, and that the 
extent of shear thinning increases with increasing solvent 
quality. However, the incremental increase in shear thin- 
ning with increasing values of z, for a given value of /3, 
appears to decrease as z increases. This can be seen 
more clearly from Figs. |3 where the asymptotic visco- 
metric ratios are plotted versus the solvent quality z, at 



FIG. 3: (a) Asymptotic non-dimensional viscosity ratio 
(?7p/?7p,o), and (b) Asymptotic non-dimensional first normal 
stress difference ratio (3 f i/\ti i o), plotted against z for two 
different values of f3. Symbols denote values obtained in the 
crossover regime, while the horizontal lines correspond to val- 
ues of the ratios in the excluded volume limit. 



two representative values of j3. In each case, the ratios 
decrease, but with decreasing rate, as the solvent quality 
increases, before appearing to level off at a final solvent 
quality independent value, in the limit of large z. In this 
excluded volume limit, consequently, universal viscomet- 
ric functions would depend only on /3, and on no other 
parameter. 

It is not possible here to explore further the approach 
of the crossover regime to the excluded volume limit, be- 
cause the present finite chain data is not sufficient to 
obtain accurate extrapolations. Obtaining such data is, 
unfortunately, currently prohibitively expensive compu- 
tationally. Nevertheless, behavior in the excluded vol- 
ume limit itself can be directly obtained, as has been 
discussed earlier. The horizontal lines in Figs. corre- 
spond to excluded volume limit values of the viscomet- 
ric ratios obtained by these means, the details of which 
will be described in greater detail shortly. It is immedi- 
ately apparent that the large z asymptote of the crossover 
regime does not coincide with the excluded volume limit. 
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Symbols are predictions of Brownian dynamics simulations, 
while lines are least square curve fits through the data. 



This is in complete contrast to the universal results, in 
steady shear flow, obtained earlier by Prakashi^ with the 
Gaussian approximation. While the general behavior in 
the crossover regime predicted by the Gaussian approx- 
imation is identical to that described above, values of 
viscometric ratios for z 3> 1 were found to smoothly 
approach those obtained in the excluded volume limit. 
This perplexing behavior, where predictions of the Gaus- 
sian approximation differ qualitatively from those of ex- 
act Brownian dynamics simulations, was also observed 
recently by Kumar and Prakash in their attempt to pre- 
dict the behavior of universal equilibrium and linear vis- 
coelastic ratios,^- Before discussing the origins of this be- 
havior in greater detail, however, it is necessary to com- 
plete the presentation of a few further results obtained 
in the excluded volume limit. 
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FIG. 5: (a) Asymptotic non-dimensional ratio ?7p/7?p,o, and 
(b) Asymptotic nondimensional ratio \E r i/\&i,o, for various 
values of /3, in the excluded volume limit. Circles represent 
Brownian dynamics simulation data, while squares represent 
results of the Gaussian approximation. The dashed lines are 
fits to the data using the Carreau-Yasuda model, Eqs. 1231 . 



B. Excluded volume limit 

Predictions in the excluded volume limit are obtained, 
as described in sec. Ill CI above, by extrapolating finite 
chain data acquired for various constant values of z*, to 
the infinite chain length limit. Figure 0] is an illustra- 
tive example, where simulations have been carried out at 
several values of N, for z* — 0.24, 0.28, and 0.35, respec- 
tively, and at a fixed value of /3 = 6, in order to acquire 
data for the ratio rj p /rj p ^. The common extrapolated 
value, in the limit N — > oo, for all the three values of z* , 
is the universal value, in the excluded volume limit, of 
Vp/Vp,0i at /3 = 6. As in the case of extrapolated values 
in the crossover regime, insufficiency in the accumulated 
finite chain data is revealed when curves for different val- 
ues of z* fail to extrapolate to a common point. In the 
present instance, this typically occurs for values of char- 
acteristic shear rate j3 ^ 10. Extrapolated data obtained 
in this manner, at various values of (3 (up to (3 ~ 10), 



for both the viscometric ratios, are displayed in Figs. 
The horizontal lines, displayed in Figs. earlier, have 
also been obtained similarly. 

Both the viscometric functions displayed in Figs. |3] ex- 
hibit a similar dependence on (3 — nearly constant values 
at small values of /3, followed by a crossover to a power 
law dependence at large values of (3. Results obtained 
earlier with the Gaussian approximation are also pre- 
sented in Figs. [S] for comparison. While the extent of 
shear thinning is under-predicted by the Gaussian ap- 
proximation for the values of j3 displayed in the figures, 
it predicts a larger slope in the power law regime, which 
implies that the shear thinning predicted by it will even- 
tually be greater than exact Brownian dynamics simula- 
tions. 

The Carreau-Yasuda model, which proposes the fol- 
lowing expressions for describing the variation of the vis- 
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TABLE I: Values of Carreau-Yasuda model parameters for 
Vp/Vpfi and ^1/^1,0- Comparison of exact Brownian dy- 
namics simulations (BDS), with the Gaussian approximation 
(GA). 





rn v 


n v 


a n 


BDS 


-0.255 ± 0.009 


1.925 ±0.100 


0.712 ±0.050 


GA 


-0.288 ±0.006 


1.675 ±0.097 


0.429 ± 0.025 








dip 


BDS 


-0.52 ±0.02 


2.395 ± 0.090 


1.312 ±0.020 


GA 


-0.750 ±0.014 


1.878 ±0.795 


0.507 ±0.012 



cometric ratios with (3, 



[i ±K 



#1,0 



[1 ± (ov, 0)"* 



has been shown previously by PrakashiS to provide a 
good fit to the results of the Gaussian approximation. 
As indicated in Figs. El with a suitable choice of the fit- 
ting coefficients a v , m n , n v etc., it also leads to a good 
fit of Brownian dynamics simulation data. The values 
of the fitting coefficients used in the two cases are dis- 
played in Table HJ The constant m represents the slope 
of the power law region at large values of (3. In line with 
the expectation by visual inspection, the value of m for 
the Gaussian approximation is larger than that for exact 
Brownian dynamics simulations. Using renormalization 
group methods, Ottinges*!^ has previously found that 
the values, m n = —0.25 and = —0.5, describe the 
power law shear rate dependence of rj p /r]pfl, and ^i/^i^, 
respectively, in the excluded volume limit. As is evident 
from the values of the power law exponents displayed in 
Table [J the predictions of Brownian dynamics simula- 
tions are in remarkable agreement with renormalization 
group results. 

The expression for Eq. 1|12|) . can be rewritten in 
the form, 



*J7,0 



(tfl/tf 



1,0J 



(Vp/Vp.o) 2 

[l±(a^ / 3)"^] (m ^ / "^ ) 
[1 + 



(24) 
(25) 



where, CAp r) ,o is the value of U^r n in the zero shear 
rate limit, and the Carreau-Yasuda model expression, 
Eq. (|23[1 . has been used for the viscometric ratios. Both 
renormalization group results, and the present Brownian 
dynamics simulations suggest that, « 2m^ w —0.5 
(see Table [Q. For /3 > 1, therefore, Eq. J2SJ) can be 
simplified to, 



S/r/,0 



(26) 



Kumar and PrakasbJi have estimated that in the ex- 
cluded volume limit, Uqj, h o = 0.771 ±0.009. Substituting 
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FIG. 6: The dependence of the universal ratio U^n on the 
characteristic shear rate j3. The open circles are predictions 
of Brownian dynamics simulations, while the squares are pre- 
dictions of the Gaussian approximation, calculated from data 
reported by Prakash. 18 The dashed curves have been obtained 
from Eq. I25H . using values for the Carreau-Yasuda model pa- 
rameters given in Table The horizontal lines represent the 
zero-shear rate predictions in the excluded volume limit, re- 
produced from ref. |3 



this result, and the values for a v and from Table^into 
Eq. I|26l) , and using standard methods for the estimation 
of error propagation, leads to = 0.568 ±0.038. Thus, 
the universal ratio Uy v decreases from a zero shear rate 
value of 0.771 ± 0.009, to a shear rate independent value 
of 0.568 ± 0.038 in the limit of large shear rates. 

Figure displays the dependence of U^, v on /3, ob- 
tained by carrying out Brownian dynamics simulations 
(open circles), and with the Gaussian approximation 
(squares). In each case, the symbols denote values that 
have been calculated using Eq. with the excluded 

volume limit values substituted for all the quantities on 
the right hand side. The dashed curves have been ob- 
tained from Eq. (|25fl , using mean values for the Carreau- 
Yasuda model parameters given in Table [I] Unlike in 
the case of exact Brownian dynamics simulations, since 
m^p 7^ 2m^ in the Gaussian approximation, U^ v does not 
reach a constant value in the limit of large [3. 

Perhaps the most celebrated result of the study of ex- 
cluded volume interactions in dilute polymer solutions is 
the observation that at large molecular weights, the equi- 
librium root-mean-square radius of gyration of the poly- 
mer coil, i? g , scales with molecular weight as a power law, 
R g ~ M u , where v is a universal exponent, independent 
of the particular polymer-solvent system. Pierleoni and 
Ryckaerti^ have shown recently, by calculating the struc- 
ture factor for bead-spring chains with up to 300 beads, 
that in the presence of excluded volume interactions, two 
distinct regimes can be identified in the structure of a lin- 
ear polymer chain undergoing shear flow. At large length 
scales, the polymer coil has a Rouse-like behavior, while 
at small length scales, where the flow field has not yet 
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distorted the isotropic coil structure, the polymer coil ex- 
hibits typical self-avoiding statistics. In the light of these 
observations we expect, therefore, that in the limit of 
infinite chain length, the scaling exponent v will remain 
unaltered from its equilibrium value, since, at any partic- 
ular shear rate, more and more polymer length scales will 
be smaller than the smallest length scale probed by the 
flow field. Recently, Kumar and Prakash^i have obtained 
the value of v at equilibrium, by exploiting a unique fea- 
ture of the solution to the excluded volume problem. We 
show here that the same procedure can also be applied 
to find the value of v in the presence of shear flow. In 
order to do so, however, it is necessary to make a few 
introductory remarks. 

While the power law scaling of R g with chain length N 
occurs only in the excluded volume limit, it is common to 
describe the dependence of R g on N, at all values of chain 
length, with an apparent power law, R g ~ 7V l ' off , with 
v e fi representing an effective exponent that approaches 
its critical value v, as N — ► oo. Schafer and cowork- 
erg 2i24 have shown that the manner in which v e g ap- 
proaches its asymptotic limit is strongly dependent on 
the magnitude of the parameter z* relative to its fixed 
point value z\. Discussions of the origin of the fixed 
point, and its significance can be found in treatises on 
renormalization group methodsAiS Basically, the exis- 
tence of power law behavior in the excluded volume limit, 
such as the one observed for i? g , is intimately connected 
to the existence of a fixed point value for the parame- 
ter z* . Schafer and coworkers have shown, using both 
renormalization group methods, and Monte Carlo sim- 
ulations, that 2/ e ff — > v on two distinct branches, (i) the 
strong-coupling branch corresponding to z* > z\ , and (ii) 
the weak-coupling branch corresponding to z* < z\ . On 
the weak-coupling branch, for increasing values of N (or 
equivalently z), v c g approaches v from below, increas- 
ing rapidly at first before approaching the asymptotic 
value very gradually. On the other hand, on the strong- 
coupling branch, i/ e g approaches v from above, decreasing 
rapidly before approaching the asymptotic value slowly. 
Kumar and Prakash2» have shown that Brownian dynam- 
ics simulations are able to capture the existence of the 
dual branched structure of the solution to the excluded 
volume problem. Furthermore, they show that this dis- 
tinctive feature of the solution implies that at fixed values 
of z and d*, a plot of v c g versus z*, for a range of val- 
ues spanning the fixed point, has a characteristic shape 
exhibiting a marked inflection point, from which one can 
obtain both the fixed point z\ , and the critical exponent 
v (see Fig. 11 in ref . l24|) . In particular, they found that 
at equilibrium, v = 0.6, and 0.28 < z\ < 0.3. 

By collecting data on R g as a function of N, at various 
constant values of the parameters z* , and f3, the value of 
f c ff, at particular values of z and d*, can be calculated 
with the help of the expression, v c q — d\nR g /d\nN. 
The symbols in Fig.0have been obtained in this manner, 
for two representative values of the characteristic shear 
rate (3. It is clear from the figure that the dual branched 
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FIG. 7: The effective exponent v e g as a function of the 
strength of excluded volume interactions z* , predicted by 
Brownian dynamics simulations, at two different vales of the 
shear rate (3. 



structure of the solution persists into the non-equilibrium 
regime, and the characteristic shape of the dependence 
of v e s on z* , observed at equilibrium, is preserved at 
finite values of the shear rate. The point of inflection, 
for the data corresponding to j3 = 0.4, clearly indicates 
that both Zf and v are unaltered from their equilibrium 
values. At = 1.75, while z\ lies in the same interval 
as at equilibrium, there is a small decrease in the value 
of v, arising perhaps due the relatively small values of 
N used to acquire the data. The results displayed in 
Fig. are consistent with the argument that, for chains 
having sufficiently large number of beads, regardless of 
the characteristic shear rate (3, there exists a local length 
scale below which the beads do not experience flow, and 
are effectively at equilibrium. 

It is appropriate now to discuss the earlier observation 
in Fig. 13 that the exact crossover functions for r] p /rj p fi 
and Wi/Wi,o, in the limit of large z, do not approach the 
asymptotic predictions obtained in the excluded volume 
limit. Kumar and PrakashSa have speculated, in connec- 
tion with a similar observation made earlier for universal 
equilibrium and linear viscoelastic ratios, that the origin 
of this behavior might lie in the fundamental difference 
between systems described by the crossover region, and 
those that correspond to the excluded volume limit. As 
pointed out earlier, since the crossover behavior has been 
obtained by keeping z constant as N — > oo (which implies 
that z* — > 0, and T — > To), points on a crossover curve 
correspond to systems, (i) that are infinitesimally close 
to the ^-temperature, and (ii) that satisfy the criteria 
for belonging to the weak-coupling branch. On the other 
hand, since the excluded volume limit behavior has been 
obtained for systems with a finite value of z* , they corre- 
spond to systems, (i) with a non-zero difference between 
T and Tg, and (ii) that belong to either the weak or the 
strong-coupling branch. Interestingly, it was shown in 
ref. |2J that, at equilibrium, the Gaussian approximation 
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does not possess the dual-branched structure of the exact 
solution. We have not examined here whether such is the 
case even at finite shear rates. Nevertheless, as pointed 
out previously, both at equilibrium and at steady state, 
the large z crossover behavior of the Gaussian approxima- 
tion coincides with the behavior in the excluded volume 
limit, unlike in the case of the exact solution. 

A number of studies aimed at predicting the orienta- 
tion angle and orientation resistance in shear flow, associ- 
ated with various tensorial quantities, as a function of the 
shear rate, have been carried out previouslyii£i22i22i2iii2ii 
The general consensus is that these quantities are univer- 
sal properties of a dilute polymer solution, in the sense 
that, at the same reduced shear rate /3, chains of differ- 
ent lengths have the same orientation angles and resis- 
tances. The universal dependence of \G an d tuq on /3, 
predicted by the present Brownian dynamics simulations 
in the excluded volume limit, will be presented shortly 
below. However, the universal shear rate dependence of 
the orientation angle Xr, and the orientation resistance 
m T (with the expected power law behavior at large val- 
ues of P), can be derived in a straightforward manner 
from the results that have been obtained above for the 
shear rate dependence of the viscometric functions in the 
excluded volume limit, namely, Eqs. 

Substituting the expressions for the viscometric func- 
tions, Eqs. 0, into the defining expressions for \r arL d 
m T , Eq. (|15|) . and using the Carreau-Yasuda model, and 
the relation between 7 and /3, leads to the following ex- 
pressions, 



cot(2xr = — ^ -, — 1 r-P 

_ U^fi [l + (a^P) n 4 m * M 
2 [1 + (a v p) n ^ m "/ n ^ 

2 (v P /Vp,o) 



U 9r ,,o (*i/*i,o) 

2 [l + K/3)"^" 1 " 7 "^ 
Uvnfl [l + ( < ty/?)"*] (m * /n * ) 



(27) 
(28) 
(29) 
(30) 



On using the relation = 2m v , Eqs. (|28|l and (|30|l sim- 
plify to the following power law dependencies of cot(2% T ) 
and m T on /3, in the limit of large 0, 



cot(2 XT ) = ^f £ (^) P 1+m » (31) 



= jf~ T P-^ (32) 

Since both renormalization group theory^S^ and the cur- 
rent Brownian dynamics simulations suggest « —0.25 
(see Table H|), Eqs. |(3TJ| and indicate that cot(2x T ) 
grows approximately as @° , while m T has a weaker de- 
pendence, growing approximately as /3 ' 25 , in the limit of 
large j3. These predictions are of course restricted to the 




FIG. 8: Universal dependence of the cotangent of the orienta- 
tion angles, (a) Xt (open circles), and (b) \G (filled circles), on 
the characteristic shear rate /3, in the excluded volume limit, 
obtained by Brownian dynamics simulations. The dashed 
curve through the Xt data has been drawn using Eq. 1281 . 
with parameter values reported in Tabled] The solid lines are 
the predictions of the Rouse model. 



influence of excluded volume interactions on these quan- 
tities, and a thorough comparison with experiment re- 
quires the incorporation of hydrodynamic interaction ef- 
fects. Bossart and Ottinger— and Cifre and de la Torre^ 
have shown previously that the presence of hydrody- 
namic interactions also leads to shear rate dependent 
orientation angles and resistances. Interestingly, how- 
ever, in contrast to the prediction obtained here in the 
presence of excluded volume interactions alone, they find 
that at large values of P, the orientation resistance m T 
decreases with increasing shear rate. As is well known, 
the presence of hydrodynamic interactions leads to the 
prediction of shear thickening at large values of /3— It 
is clear from Eq. (|32|l . that a positive value of the expo- 
nent rrifj, would explain the observation of Bossart and 
Ottinger, and Cifre and de la Torre. 

Figures El display the shear rate dependence of the 
orientation angles obtained here by Brownian dynamics 
simulations in the excluded volume limit. In the case 
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FIG. 9: Universal dependence of the orientation resistances 
m T (squares) and rac (circles), on the characteristic shear 
rate f3, in the excluded volume limit. The line through the m T 
data has been drawn using Eq. 13UH . with parameter values 
reported in Table Q The horizontal lines are the predictions 
of the Rouse model. 



of cot(2x T ), values for rfp/r/po and ^1/^1,0, obtained at 
various values of (3, in the excluded volume limit, have 
been substituted on the right hand side of Eqs. J37J, with 
U* n ,o = 0.771 ± 0.009. On the other hand, cot(2xc) has 
been evaluated by accumulating finite chain data using 
Eq. (|16fl . at a constant value of z* = 0.28 (the fixed point 
value), followed by extrapolation to the long chain limit. 
Comparison with the Rouse model (which predicts a lin- 
ear dependence on (3 in each case) , clearly indicates that 
the presence of excluded volume interactions leads to a 
weakening of the alignment of the polymer coil in the 
flow direction. This has been observed previously both 
theoretically 15 and experimentally. 38 Since, at the same 
shear rate /3, cot(2xc) > cot(2x T ), the tensor G is more 
easily oriented than t p . 

Figure displays the universal dependence of the ori- 
entation resistances mg and m T , in the excluded volume 
limit, on [3. As in the case of Xn the orientation resis- 
tance m T has been obtained by substituting known values 
on the right hand side of Eq. (|29(l . On the other hand, 
ma has been obtained from the Brownian dynamics sim- 
ulation data for cot(2x(j), displayed in Fig-Efb), by using 
the defining expression, vtlq — f3/ cot(2xG)- It is evident 
from the figure that both ma and m T increase with /?, 
indicating that it gets increasingly difficult to orient the 
polymer coil with increasing shear rate. This has been 
previously predicted^ and observed experimentally at 
relatively low values of [3~L The fact that G is more easily 
oriented than r p can also be seen from the behavior of the 
orientation resistances rriQ and m T , since m r > mg, for 
all values of (3. As pointed out previously by Bossart and 
Ottinger^ flow birefringence experiments (which mea- 
sure Xr)i reflect orientational ordering on length scales 
smaller than those probed by light scattering experiments 
(which measure xg)- Since smaller length scales are less 



affected by flow, they exhibit greater resistance to orien- 
tational ordering, leading to the observed relative mag- 
nitudes of me and m T . 

It is of particular interest to examine the universal val- 
ues of the orientation resistances in the limit of zero shear 
rate, namely, T71g,o and m T; o, since a number of studies 
have been confined to the low shear rate regime. Bossart 
and Ottinger, 32 and Pierleoni and RyckaerfeS have ar- 
gued previously that the effect of excluded volume inter- 
actions on the orientation resistance is negligible. Bossart 
and Ottinger^i have found, from careful experiments, 
that in a good solvent (polystyrene in bromobenzene) , 
«V,o = 3.74 ± 0.28, compared to the value of 3.13 ± 0.42 
in a theta solvent (polystyrene in 4-bromo-a-benzyl alco- 
hol) . Given the magnitude of the error bars, it can only 
be concluded that m Tj o is probably larger in a good sol- 
vent than in a theta solvent, ft must be borne in mind 
that both these values reflect the influence of hydrody- 
namic interactions. In the present instance, where only 
the influence of excluded volume interactions is consid- 
ered, the value of m Ti o can be obtained from the following 
expression, 

m r,o = -rr— (33) 

which is derivable from Eq. I|29|l in the limit of zero shear 
rate. Using the previously established exact Brownian 
dynamics simulations value of U^, n .Q = 0.771 ± 0.009, 24 
leads to, m T ^ — 2.594 ±0.030, which is marginally larger 
than the Rouse value of 2.5. Using non-equilibrium 
molecular dynamics simulations of bead-spring chains 
with FENE springs, Aust et al.^i have estimated m r ,o = 
2.4 ± 0.3. Interestingly, in spite of the fact that the 
solvent is treated explicitly in these simulations, and 
consequently hydrodynamic interactions are taken into 
account, the predicted value is not significantly differ- 
ent from the pure excluded volume prediction obtained 
here. It is worth noting that Bossart and Ottinger, 32 
who used the Gaussian approximation, and de la Torre 
and co-workersj22i^ who used Brownian dynamics sim- 
ulations with hydrodynamic interactions, report the val- 
ues: 777 T ,o ~ 3.57, and m Ti o 3.5, respectively — both 
of which are substantially different from the prediction 
of the Rouse model. While an improvement in solvent 
quality seems to increase 777 Ti o fractionally, our results 
seem to suggest that it decreases 777^0 fractionally. By 
fitting the ma data in Fig. [5] with a least squares third 
order polynomial, and extrapolating to (3 — > 0, we find 
that in the excluded volume limit, mc,o = 1-69 ± 0.03, 
which is slightly smaller than the Rouse prediction of 
wiQfi = 1.75. The value obtained here is in good agree- 
ment with the explicit solvent non-equilibrium molecu- 
lar dynamics simulations of Aust et al.— and the results 
of renormalization group calculations, 32 both of which 
predict mc,o — 1-7. Since these studies account for hy- 
drodynamic interactions, it is clear that hydrodynamic 
interactions have no effect on mc,Q- Pierleoni and Ryck- 
aert have also observed a similar pattern when comparing 
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their explicit solvent molecular dynamics simulations 39 
with their Brownian dynamics simulations (in which only 
excluded volume interactions were taken into account) ^ 
It is well known that hydrodynamic interactions have a 
significant influence on dynamic properties, but have no 
influence on static properties, since the equilibrium dis- 
tribution function is unaltered in their presence. This 
suggests, perhaps not surprisingly, that mc.o is static 
property, while to t .o is a dynamic property. 

IV. CONCLUSIONS 

Results of a detailed Brownian dynamics simulations 
study, of a polymer solution undergoing steady shear 
flow, have been presented. The polymer molecule has 
been modelled by a bead-spring chain, and a narrow 
Gaussian excluded volume potential, that acts between 
pairs of beads, has been used to mimic the influence of 
solvent quality. 

Material properties have been shown to become in- 
dependent of the range of excluded volume interactions 
d*, and the number of beads N, in the limit of large 
N. Furthermore, it has been found that master plots 
are obtained when, (i) the influence of the strength 
of excluded volume interactions is interpreted in terms 
of the solvent quality parameter z = z*vN, and (ii) 
the dependence on the shear rate 7 is interpreted in 
terms of the characteristic non-dimensional shear rate, 
/3. In this work, we have explored the universal depen- 
dence of the non-dimensional viscosity ratio (?7 P /?7p, o)j 
the non-dimensional first normal stress difference ratio 
(^i/^i^o), the orientation angles \t and \G, an d the 
orientation resistances m T and tog, on the characteristic 
shear rate /3 and the solvent quality z. 

The extent of shear thinning has been found to in- 
crease as z increases (Figs. |2J. However, the incremen- 
tal increase in shear thinning saturates as z increases 
(Figs. |3J). The existence of universal viscometric ratio 
versus shear rate curves, independent of all model pa- 
rameters, in the excluded volume limit, z — > 00, has 
been verified (Figs. EJ). The shear rate dependence of 
both the universal viscometric ratios has been found to 
be well described by the Carreau-Yasuda model, with 
power law decay exponents at large values of j3, that are 



in close agreement with values predicted by renormaliza- 
tion group theory (Table. 0}. 

The asymptotic values of the exact crossover functions 
for (r? p /??p,o) and (\f r i/\f'i,o), m the limit of large z, are 
found not to approach the universal predictions obtained 
in the excluded volume limit, unlike previous predictions 
of the Gaussian approximation (Figs.[3J). 

The shear rate dependence of the universal ratio U^^ 
has been obtained in the excluded volume limit. The 
ratio is found to decrease from the zero shear rate value 
of 0.771 ± 0.009, to a constant value of 0.568 ± 0.038, in 
the limit of large shear rates (Fig. . 

The dual branched structure of the solution, eluci- 
dated recently by Scafer and co-workers and shown 
to be captured by Brownian dynamics simulations at 
equilibrium^ has been found to persist into the non- 
equilibrium regime. As in the equilibrium case, this dis- 
tinctive structure has been exploited here to obtain an es- 
timate of both the fixed point of the strength of excluded 
volume interactions, and the critical exponent v. Both 
quantities are found to remain unaltered in the presence 
of shear flow (Fig.0). 

The universal dependence of % T , XGi m T and tog, on 
(3, in the excluded volume limit, has been obtained. The 
presence of excluded volume interactions has been shown 
to weaken the alignment of polymer coils in the flow di- 
rection (Figs. (SJ. Of the two tensors, the radius of gy- 
ration tensor G is found to be more easily oriented than 
t p (Fig. El). The power law shear rate dependence of 
Xr (~ /3°' 75 ) and m T (~ /3°' 25 ), in the limit of large 
shear rate, has been obtained by exploiting the Carreau- 
Yasuda model fit of the universal viscometric functions. 
The presence of excluded volume interactions leads to 
a marginal increase in to Tj o, while marginally reducing 
m G ,o- 
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